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Since the prediction of a new topological state of matter in graphene, materials acting as topolog- 
ical insulators have attracted wide attention. Shortly after the theoretical proposal for a mercury 
telluride (HgTe)-based two-dimensional topological insulator, the observation of the quantum spin 
Hall effect and non-local edge transport brought compelling experimental evidence for quantized 
conductance due to edges states. The spin orientation and propagation direction of such helical 
edge states are inherently connected, providing protection against backscattering. However, these 
features conversely render controlled spin operations such as spin switching difficult. Here we there- 
fore propose constrictions as connectors between opposite edge (and spin) states in HgTe. We 
demonstrate how the coupling between edge states, which overlap in the constriction, can be em- 
ployed both for steering the charge flow into different edge modes and for controlled spin switching. 
This gives rise to a three- state charge and spin transistor function. Unlike in a conventional spin 
transistor, the switching does not rely on a tunable Rashba spin-orbit interaction, but on the energy 
dependence of the edge state wavef unctions. Based on this mechanism, and supported by numerical 
transport calculations, we present two different ways to control spin- and charge- currents depending 
on the gating of the constriction. 



Realizations of topological insulators-^ with their ex- 
ceptional electronic and transport properties exist both 
in three dimensions, e.g. the surface states of BiSe^ or 
Bi2Te3 4,5 , and in two dimensions with HgTe as promi- 
nent example. The topological properties of HgTe/CdTe 
heterostructures are attributed to the combination of an 
inverted band structure of the HgTe layer (the light- 
hole band as conduction band and the heavy-hole as the 
topmost valence band) and the regular band structure 
of CdTe. Quantized transport along the HgTe bound- 
aries can be conveniently explained by an edge channel 
pictured Two states with opposite spin orientation prop- 
agate along opposite device edges in the same direction 
and thus lead to a quantized conductance of 2e 2 /h. Due 
to the spatial separation of the spin-states the spin-orbit 
coupling is suspended, and the system geometry can be 
employed for spin selection^. 

Spin-selectivity is also a crucial element of the Datta- 
Das spin transistor proposal^, where charge flow is con- 
trolled electrically through the gate-dependent Rashba 
spin orbit interaction^ (SOI) in a conventional two- 
dimensional semiconductor heterostructure placed in be- 
tween ferromagnetic contacts. Its realization, however, 
turns out to be difficult owing to spin relaxation in 
the semiconductor heterostructure and interfacial effects 
such as the conductivity mismatch between the different 
materials 10 . HgTe-based topological insulators appear to 
be promising candidates for spin processing devices since 
they also can be gated and exhibit considerable SOI but, 
on the contrary, are composed of a single material class 
only. Moreover, the one-dimensional (Id) nature of their 
edge states suppresses orbital effects present in bulk con- 
ductors, leading to high spin polarizations and to a much 
better (spin) switching quality. 

To our knowledge there have been only a few proposals 
for spin-transistors based on two-dimensional topologi- 
cal insulators. Two of them rely on spin switching with 



a magnetic field at a pn-j unction^ or in an Aharonov- 
Bohm interferometer—. Recently it has further been 
suggested, within a phenomenological model, that sep- 
arate gating of the two branches of an Aharonov-Bohm 
interferometer allows for manipulating charge- and spin- 
transport^. By contrast, our present proposal relies on 
an electrical operation by gates on a single HgTe nano- 
constriction, which up to now was only considered for 
charge current switching^. 

In this manuscript we demonstrate how topologi- 
cal edge states can be selectively switched in an elon- 
gated constriction etched out of a HgTe heterostructure, 
leading to an integrated three-state charge- and spin- 
transistor of high fidelity. An incoming spin-polarized 
(upper edge) state can either be reflected back to the 
lower edge, as shown in Fig. QJa), or transmitted through 
the narrow part. Back-reflection into the opposite spin 
channel at the same edge is forbidden, as the absence of a 
magnetic field implies time-reversal symmetry^. Within 
the constriction the SOI between the edge channels is re- 
activated due to a finite overlap between right moving 
edge-states on upper and lower side, giving rise to spin 
precession that is tunable by a gate. This allows for steer- 
ing the spin orientation of the electrons which leave the 
constriction, and thereby their further path: An incom- 
ing spin-up state will leave the system either by swapping 
the edge (with a simultaneous spin- flip), as shown in Fig. 
Gib), or by remaining in its spin and edge state as shown 
in Fig. Die). 

In the following, we first introduce a Id model which 
explains all relevant spin and charge transport features 
arising from the constriction. Subsequently, we check the 
predictions of this model by numerical quantum trans- 
port calculations which include random impurity poten- 
tials and rough confining potentials to account for ma- 
terial inhomogeneities and etching imperfections. Fur- 
thermore, we analyse the peculiar switching properties 
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by means of a top-gate and side-gates acting on the con- 
striction. 



We describe the electronic properties of the underlying 



HgTe heterostructure by the four band HamiltoniarAi^ 
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FIG. 1. (color online) Local density of states and band 
structure of a HgTe constriction for different switching con- 
ditions - a-c, Spin-resolved local density of states for charge 
carriers entering a disordered constriction from the upper left 
edge. The colour code indicates the spin polarization of the 
state, with blue representing the upper subblock (t) and or- 
ange the lower subblock (1). a, The chemical potential in the 
narrow region is chosen within the confinement induced gap 
([i = 7 meV), leading to a perfect reflection of the incoming 
state, b, For an energy closely below the gap (ji = 4 meV) the 
effective spin-orbit coupling is pronounced enough to produce 
a precession from spin-up to spin-down and the associated 
switching from the upper (left) to the lower (right) edge, c, 
At low energies, distinctly below the gap (/i = —5 meV), 
the edge states are predominantly localized at the respective 
boundary in the constriction, leading to a reduced effective 
SOI, and hence the incoming spin-up states leave the con- 
striction unrotated. d, Sketch of the HgTe band structure for 
the sequence of a wide (left panel, W = 1000 nm), narrow 
(middle panel, W = 100 nm) and wide (right) lateral confine- 
ment. Dashed lines mark the energies of the local density of 
states in panels a-c. 
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where k± = k x ± ik y and k = k 2 + k 2 . This Hamiltonian 
contains the commonly used time-inverted 2x2 blocks 
for the composite states of the heavy-hole and electron 
bands 16 . Additionally, we take into account the leading 
order SOI terms A and R, due to bulk-inversion asym- 
metry (BIA) and structure-inversion asymmetry (SIA)i£, 
respectively. (See the supplementary material for the val- 
ues used for A, R and the further material parameters 
A, B, D and M.) 

Based on the Hamiltonian (pQ) we first derive an effec- 
tive Hamiltonian for an infinite strip of constant width 
W. We chose the lead to point in ^-direction and search 
for the transversal eigenfunctions ip(y) which separately 
fulfill the boundary conditions ip(y < 0) = at the lower 
side and ijj(y > W) = at the upper side. For a very 
wide confinement one can neglect the influence of the 
opposite boundary on the edge states, leading to a full 
spin-polarization. Then the resulting states can be clas- 
sified by their subblock into up (t) and down (^), as well 
as their propagation direction into right movers (+) and 
left movers (— ). Accordingly, the right moving state of 
the upper subblock is given by 

1>+(y) oc ( e *i(y- w ) -e x ^y- w >) (l, -£+, 0, 0) (2) 

and the right moving state of the lower subblock by 

tl>+{y) oc (e" A ^ - e~ x iy) (0, 0, 1, £ + ) , (3) 

with two different decay exponents , > and £ + 
the weight of the second spinor entry in the respective 
subblock. Then the left moving states are obtained by 
complex conjugation and by substituting y — ^ W — y and 
k x — >• — k x which corresponds to inter-changing A^ 2 — >• 

\~l 2 and £ + — >> £~ . (The functions are described in detail 
in the supplementary material.) From these properties 
we can derive an effective Id Dirac Hamiltonian H f 
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c =p aa x k x for a sing le spin subblock with a velocity a = 
A y/(B 2 -D 2 )/B 2 and an energy offset c = -DM/B. 
For a wide strip this is in perfect agreement with the full 
band structure in the vicinity of the band crossing shown 
in the left and right panel of Fig. QJd) for W = 1000 nm. 

For decreasing width W, the edge states at opposite 
boundaries start to overlap, leading to a mass like gap 
in the Id Hamiltonian. By invoking simultaneously the 
boundary conditions for the upper and lower side the size 
of the effective mass gap is given by^ 
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Additionally, the suppression of the SOI for the edge 
states is suspended for small W. Neglecting the rapidly 
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decaying terms proportional to e x ^ y in the wavefunction, 
the overlap due to BIA can be stated as 
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The effect of SIA on the edge states within the band gap 
is negligible small. Combining the effective mass m of 
Eq. (|H) and the effective SOI S ± of Eq. ©, we can 
compose a Id effective Hamiltonian 




(6) 



with 5p/ m = (£ + ± S~)/2. The band structure of i^ e ff 
exhibits a mass gap and a variable, energy dependent 
effective spin-orbit splitting which is strongest close to 
the avoided band crossing shown in the middle panel of 
Fig. DJd) for a constriction of width W = 100 nm. Note 
that the SOI of the model and the corresponding splitting 
in the band structure is slightly overestimated compared 
to the result of the full Hamiltonian, due to the restriction 
to only four basis functions. 

In the following, we analyze within this model the 
transport properties for a constriction interconnecting 
two bulk-like regions in an H-shaped HgTe heterostruc- 
ture, as depicted in Fig. Efa). Within the bulk bandgap 
the transport is exclusively carried by edge states. Ac- 
cordingly, there are three different paths entering at lead 
A and leaving the system at lead B (green), C (orange) 
or D (blue). 

The Id model for this setting equals a Dirac equation 
with a position-dependent mass potential, if we neglect 
SOI. For an abrupt change in width, this model can be 
solved analytically 19 leading to a pronounced drop of the 
transmission from A to D in the energy region of the con- 
finement induced gap as well as strong Fabry-Perot like 
transmission resonances for energies below the gap. In 
our Id model calculations we use a more realistic smooth 
variation of the width given by two Fermi functions, 



W(x) = W 



Wp -W Wp-W \ 

l + e (x-L/2)/5L ~ l + e {x+L/2)/5L J ' 

(7) 

Here W is the width of the constriction of length L, Wo is 
the width outside the constriction, which is chosen wide 
enough to ensure a gapless Dirac spectrum, and SL de- 
notes the length of the transition from Wo to W, see 
Fig. [2^a). The transmission through the smooth con- 
striction, as shown in Fig. Efb), is calculated using a 
discretized transfer-matrix approach. As a result the 
Fabry-Perot type resonances are strongly reduced while 
the transmission gap remains unaffected. This simple 
model can readily be used to describe switching 14 of the 
output current between lead B or D upon tuning the 
Fermi energy, as shown in Fig. [2fb). 
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FIG. 2. (color online) Quantum transmission through a con- 
striction: one-dimensional model vs. full quantum mechan- 
ics - a, Sketch of the constriction geometry. The coloured 
paths illustrate the differently scattered edge channels enter- 
ing through lead A and leaving the system at lead B, C or 
D. b, Transmission through a constriction neglecting SOI cal- 
culated for the Id model Hamiltonian, see text. Different 
colours mark transmission paths according to colour coding 
in (a), c, Energy dependence of the effective SOI, Eq. J5}. 
d, Transmission through a constriction with SOI calculated 
by means of the Id model Hamiltonian ©. e, f, Transmis- 
sion obtained by a full numerical wave-packet calculation for 
(e) a perfect constriction and (f) a realistic constriction with 
random impurity potential (Uq = 2meV) and rough walls 
(W r = 20nm), the same system as used in the calculation 
shown in Fig. [TJa-c). Geometric properties of the constric- 
tion: width W 100 nm, length L 1900 nm (SL 50 nm), 
width of bulk parts Wo — 1000 nm. 



However, neglecting the SOI forecloses the possibility 
of a unique spin-flip mechanism as shown in the follow- 
ing. With SOI the right moving spin-up and spin-down 
states hybridize within the constriction to a symmetric 
and an antisymmetric composite state with a difference 
in the wave vectors of Afc ~ \5 + \/a. Since the overlap 
strongly depends on the extent of the edge-states, which 
again is mainly determined by the decay exponents Xf , 
the effective spin-splitting energy [in Eq. (j5j)], has 
a pronounced energy dependence as shown in Fig. EJc). 
For low energies (larger Xf) both states are very closely 
bound to opposite edges. Thus the effective SOI is very 
weak and the states preserve their initial spin while trav- 
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eling through the constriction [see also Fig.QJc)], leading 
to a perfect transmission from lead A to D as shown by 
the blue curve in Fig. [2(d). At higher energies, close 
to the confinement induced bandgap, the spin splitting 
is pronounced and already visible in the band structure 
shown in the center panel of Fig. QJd). Hence the states 
undergo a spin precession when traversing the constric- 
tion. Upon tuning the energy, i.e. |J ± (E')| [see Fig. [1(c)], 
the precession frequency can be steered. This spin rota- 
tion leads to a finite transmission from lead A to C, shown 
as the orange line in Fig. [2(d). For a given constriction 
length there is an energy (around 2 meV) where the in- 
coming spin is flipped leading to nearly perfect transmis- 
sion into lead C. Note that the three different transmis- 
sion probabilities displayed in Fig. [2(d) show well sepa- 
rated and pronounced maxima up to unity for the various 
paths. As a result the Id model suggests that a HgTe con- 
striction is perfectly suited to function as a three state 
spin-orbit transistor, switching with excellent on-off ra- 
tio between the outgoing leads B, C and D. At the same 
time, this system allows for controlled spin swapping, i.e. 
when choosing path A to C. 

So far we employed aid model to demonstrate the vari- 
ous charge and spin switching properties of a HgTe-based 
nano-constriction. In the following we study the robust- 
ness of these effects for a realistic setting. To this end 
we calculate the conductance governed by the four band 
Hamiltonian (pQ) through constrictions in two-dimensions 
with additional random impurity potentials and rough 
walls. We have extended an efficient numerical method 
to calculate electronic transport by means of the time- 
evolution of wave-packets 20 to arbitrary spin-orbit cou- 
pled systems. The propagation is calculated by means 
of an expansion of the time-evolution operator in Cheby- 
chev polynomials 2 ^, ensuring negligible numerical errors 
for long propagation times 22 . A detailed description of 
this algorithm and the considered models for impurity 
potential and edge roughness can be found in the sup- 
porting material. 

To check the validity of the Id model we first apply 
the wave-packet transport algorithm to a clean system 
with the same parameters. The resulting energy depen- 
dent transmissions into leads B, C and D are shown in 
Fig. [2(e). As mentioned before, the Id model Hamilto- 
nian [Fig. [2(d)] overestimates the effective SOI. Conse- 
quently the spin-precession is reduced in the full calcu- 
lation, leading to maximally a bit more than half a spin 
precession for the chosen geometry. Note that, as in the 
Id model, there is a distinct energy at which the spin pre- 
cession produces an exact spin- flip, leading to a perfect 
transmission from lead A to lead C [shown as the orange 
line in Fig. [2(e)]. Also the transmission from lead A to 
B (green line) in the energy region of the confinement 
induced mass gap and the transmission to lead D (blue 
line) is in good correspondence with the model. 

In order to approach realistic, experimentally avail- 
able conditions, we add an impurity potential with Uo = 
2meV and a wall roughness of W r = 20 nm and cal- 



culate the transport through such a system. The re- 
sulting energy dependent transmissions are summarized 
in Fig. [2(f) . Due to impurity induced energy variations 
within the constriction and fluctuations in its width, the 
mass gap is enhanced as visible in the transmission from 
lead A to B (green line). Most notably, the strength of 
the spin-flip mechanism is maintained, although the effi- 
ciency of the spin flip process is slightly reduced by the 
impurity potential and the wall roughness: The spin flip 
transmission [from lead A to C, orange line in Fig. [2(f)] 
no longer reaches T = 1 [see also Fig. [1(b)]. Nevertheless, 
this reduction only amounts to 20% for a very strong per- 
turbation as used in these calculations. Consequently we 
conclude that the switching properties of such a device 
are robust against electrostatic impurities and persist in 
non perfectly etched heterostructures. 

In the following we consider the possibility of con- 
trolled switching between edge currents by an additional 
gate. We model local gating that has been proven ex- 
perimentally feasible 7 , by a position dependent potential 
which is switched on outside of the confined region and 
spans a length of 3000 nm [including the constriction, 
see Fig. [3(a)]. Furthermore, a random impurity poten- 
tial as well as edge roughness are again considered, as 
specified above. By means of our wave-packet algorithm 
we calculate the quantum transport through the device 
as a function of both Fermi energy and gate voltage. 
The resulting transmissions into leads B to D are shown 
in Fig. [3(b-d), where large (small) transmissions are de- 
picted by dark (bright) colours. Note that all plots show 
universal conductance fluctuations for energies outside 
the bulk bandgap [marked by dashed lines in Fig. [3(b- 
d) and given approximately by Ep — V g < —\M\ and 
E F + V g > \M\]. Within the bulk bandgap the device 
exhibits the switching properties. The transmission from 
lead A to lead B, which is close to 1 in the confinement 
induced mass gap, is shown in Fig. [3(b); the transmission 
into lead C, accompanied by a spin-flip, in Fig. [3(c); and 
the transmission to lead D without spin- flip in Fig. [3(d). 
The Fermi energy within the whole system is pinned to a 
certain value and can be globally tuned e.g. by means of 
another back-gate. Assuming for example Ep = 0, the 
transmission can be steered between the three different 
leads by changing the gate voltage locally on top of the 
constriction and thereby changing the effective SOI. In 
this case a perfect transmission to lead B is achieved for 
V g = — 7meV. For V g = — 4meV the transmission to lead 
C is maximal, whereas for V g = 5meV the transmission 
to lead D approaches unity. The fact that the spin flip 
process is slightly reduced [see Fig. [3(c)] mainly stems 
from an average potential difference between the upper 
and the lower wall, which can be attributed to the elec- 
trostatic impurity potential. This leads to a deviation 
between the momenta of the right moving states from 
the two different subblocks, and an associated reduction 
of the transition (and spin-flip) probability. 

A similar effect occurs when we use two side-gates close 
to the constriction with an opposite applied voltage V s , 
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E F (meV) E F (meV) 

FIG. 3. (color online) Top-gate- and side-gate-switching of a HgTe constriction - a, Scheme of the constriction (width 
W — 100 nm, length L = 1900 nm) coupled to bulk parts (width Wo — 1000 nm) and subject to a top-gate (light red area). For 
the transport calculation a random potential of Uo = 2 meV and a wall roughness of W r = 20 nm are additionally included. The 
impurity potential and wall roughness are visualized in Figs.[TJa-c). Panels b-d show the (colour coded) quantum transmission 
from lead A to the lead B (b), C (c) and D (d) as a function of the Fermi energy Ef and the top-gate voltage V g for the 
constriction as depicted in panel a. e, Scheme of a smaller constriction (width W = 60 nm, wall roughness W r = 10 nm, length 
L = 900 nm coupled to bulk leads of Wo = 1000 nm) with two side-gates, e-f, Transmission from lead A to lead B (f), C (g) 
and D (h) as a function of Ef and side gate voltage V s . Panels (b-d) and (f-h) feature a spin-transistor switching. 
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as sketched in Fig. [3fe). This leads to different chemi- 
cal potentials for the two spin channels close to the up- 
per and lower side of the constriction. The results of a 
corresponding conductance calculation are summarized 
in Fig. [3jf-h). Due to the different gating, the mass 
gap now shows up as a vertical stripe (Ep > 5 meV), 
and the constriction is isolating, independent of the gate 
voltage T/g, as shown in Fig. [3ff). If the constriction is 
tuned into the state conducting charge from left to the 
right (Ep < 5 meV), it works as a spin transistor con- 
trolled through the side-gate voltage V s : For low |V^| the 
states entering the device undergo a spin-flip within the 
constriction whereas for larger |V^| the momenta of dif- 
ferent spin states differ sufficiently to suppress the spin 
precession. For an illustration, compare the transmis- 
sion from lead A to C in panel (g) with the transmission 
to lead D in panel (h). As before, the calculations are 
performed under realistic conditions, including a random 
impurity potential and wall roughness which reduce the 
spin flip probability. Nevertheless, the remaining signal is 
still pronounced enough to be observed experimentally in 
a single configuration. The calculations (f-h) were per- 
formed for a narrower constriction (W = 60 nm) and 
the same amount of disorder to demonstrate that the 
switching functionality is robust against down scaling to 



a regime of a few 10 nanometres. In view of Eq. (J5j), 
the effective SOI increases with decreasing width, allow- 
ing for faster spin precession and shorter constrictions 
[900 nm in Fig. EJe-h)]. 

We have shown that a constriction joining together 
edge channels of the topological insulator HgTe acts as a 
transistor with unique charge and spin switching proper- 
ties. These are robust against disorder and edge rough- 
ness as present in experiments. Mediated through an 
effective spin orbit coupling arising in the constriction, a 
local top gate enables switching between the edge states, 
while side gates allow for pure spin transistor action. 
Such nano-constrictions may serve as building blocks and 
connectors for more complex spin- and charge-selective 
edge channel networks based on topological insulators. 
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